Thermodynamic Properties of Fermi Systems with Flat Single-Particle Spectra 
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The behavior of strongly correlated Fermi systems is investigated beyond the onset of a phase 
transition where the single-particle spectrum £(p) becomes flat. The Landau-Migdal quasiparticle 
picture is shown to remain applicable on the ordered side of this transition. Nevertheless, low- 
temperature properties evaluated within this picture show profound changes relative to results of 
Landau theory, as a direct consequence of the flattening of £(p). Stability conditions for this class of 
systems are examined, and the nature of antiferromagnetic quantum phase transitions is elucidated. 
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Low-temperature properties of strongly correlated 
Fermi systems exhibit features not inherent in Landau 

Fermi liquids QSHSEEHHHIIS ° f Sp °" 
cial importance is the scaling behavior x~ (2\ H) = 
X _1 (0, H) + T a F(H/T) of the inverse magnetic suscep- 
tibility, observed so far only in heavy fermion metals 
|H [E E3 El in weak, sometimes tiny external mag- 
netic fields H. This feature rules out the collective, spin- 
fluctuation scenario advanced in Refs. fill If3 to explain 
the non- Fermi-liquid (NFL) behavior of these systems, 
while providing evidence for the direct relevance of single- 
particle degrees. Earlier work within this picture |l3l| has 
evaluated thermodynamic properties on the "metallic" 
side of a phase transition associated with the rearrange- 
ment of the Landau state in strongly correlated Fermi 
systems. Here we focus attention on the situation be- 
yond the critical point where the Landau state becomes 
unstable. 

This state is known (e.g. from Refs. H3. ll5|) to lose its 
stability at a density n = n& where a bifurcation point 
p = Pb emerges in the equation 



£(p,n,T = 0) = 



(1) 



which ordinarily has only the single root p — pf- (Here 
£(p) = e{p) — n is the single-particle energy measured 
from the chemical potential /i.) We shall focus on the 
case \pf — Pb\ ^ Pf, where the critical single-particle 
spectrum has the form 



£{p,n b ,T = 0) ~ (p-Pb) 2 (p-PF) 



(2) 



Thus if pb coincides with pp, the effective mass M*(n) 
diverges at rit,; in the general case pb ^ pf, the Lan- 
dau state loses its stability before M* becomes infinite. 
Beyond the critical density n&, Eq. (JJJ possesses two ad- 
ditional roots pi and p2 with p\ < pb < pi- The spectrum 
£,fl(p,T — 0), evaluated with the FL momentum distri- 
bution hfl(p) = 0{pf — p), then takes the form 



If pb 7^ pp, the roots p\, P2 are both located either in 
the interior of the Fermi sphere or both outside it. If 
Pb = Pf, then pi < pp < P2- In all these cases, the 
occupation numbers ufl(p) are rearranged. As a rule, 
the Fermi surface becomes multi-connected, but at T — 
0, the quasiparticle occupation numbers continue to take 
values or 1. Hence the Landau-Migdal quasiparticle 
picture holds, with n(£) = 1 for £ < and otherwise. 
Consider first the case p\ < P2 < Pf- Then according 
to Eq. (J2J), the single-particle states remain filled in the 
intervals p < p\ and P2 < P < Pf, while the states with 
Pi < p < P2 are empty. We call this new phase the 
bubble phase. If pb = Pf, then p\<pp and P2>Pf, and 
the states with p<p\ and with pf<P<P2 are occupied, 
while those for pi<p<pp are empty. Again one deals 
with the bubble phase. 

At this point, we observe that the solution © is 
not self-consistent, since the spectrum is evaluated with 
tifl{p) while the true Fermi surface is doubly-connected. 
Following Ref. Q, we consider the feedback of the rear- 
rangement of ufl{p) on the spectrum £(p) in the bubble 
phase based on the Landau relation |16l Il7| 



de{p) 



_P_ 

M 



/(P,Pi) 



dn(pi) 



(4) 



dp M J '^"-^ dpi 

where / is the scalar part of the Landau interaction func- 
tion and n(p) — [1 + exp(£/T)] -1 is the quasiparticle 
momentum distribution. For the electron liquid within a 
solid, p/M is to be replaced by de^/dp, where e° is the 
LDA electron spectrum. To date, Eq. J2J has been solved 
only in 3D Fermi systems with functions / depending on 
q = \p — pi | . Despite the diversity of forms assumed for 
f(q), the resulting spectra and momentum distributions 
bear a close family resemblance. This robustness is illus- 
trated in Figs, n and [21 which display results 0] from 
solution of Eq. J3J for the interaction function 

f(q) = Xl/[(q/pp) 2 + f3 2 1 } , (5) 
and in Figs. Eland 01 which present results for 



£,fl(p, n,T = 0) ~ {p - Pi){p - P2){p - Pf) 



(3) 



f{q) = \2/[{{q/2p F ?-l? + f3l 



(0) 
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How do the bubble solutions of Eq. J2J evolve un- 
der variation of T? At extremely low T < Tfl ~ 
(P2 — Pi) 2 /M, these states are described by FL theory 
with the enhanced effective mass M* ~ Mpp/(p2 — Pi)- 
Heating above Tfl results in their dissolution (see Figs.Q] 
and |2J). With further increase of T, the spectrum £(p) 
becomes smoother, and in the region of a new critical 
temperature Tz, a flat portion £ ~ appears over an 
interval [pi,pf] surrounding the point pf, as shown in 
the left panels of Figs. 03 and The presence of this flat 
portion of is a si gna ture of the phenomenon called 
fermion condensation |l8l ll9L l2fj| . Since — e(p) — fx 
and e(p) = SEq/ '6n(p), the equality £ = can be rewrit- 
ten as a variational condition 



SEq 

5n(p) 



Pi <P<Pf , 



(7) 



with^o = Y^ P £ p n (P) + lE P , Pl /(p-Pi)«(p)«(Pi) an d 
e p = p 2 /2M. The solution no(p) of Eq. Q, or equiva- 
lently of Eq. J7J), is a continuous function of p with a 
nonzero derivative dno / dp. The set of states with £ (p) — 
is called the fermion condensate (FC), since the asso- 
ciated density of states contains a Bose-liquid-like term 
•qn5{e). The dimensionless constant rj ~ (pf —pi)/pp is 
identified as a characteristic parameter of the FC phase. 

It has been shown [2(| that the FC "plateau" in 
has a small slope, evaluated by inserting no(p) into the 
above Fermi-Dirac formula for to yield 



i;(P,T>T z )=T\ii 



1 - nn(p) 
na(p) 



Pi <P <Pf ■ (8) 



As indicated in the bottom-right panels of Figs. 03 and 0] 
at T > T z the ratio £(p)/T is indeed a T-independent 
function of p in the FC region. The width £(p/) — = 
£/— £i of the FC "band" appears to be of order T, almost 
independently of 77 > 7y m i n ~ 1CP 2 . Thus at r\ > ?7 m i n the 
FC group velocity is estimated as 



V d P )t VPf 



Pi <p < pf 



(9) 
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FIG. 1: Single-particle spectrum f(p) in units of 10~ 3 £%, 
where e% = p 2 F /2M (top panel), occupation numbers n(p) 
(middle panel), and d^/dp in units of 10~ s v ( f, where v% = 
Pf/M (bottom panel), plotted versus p/pf at four color- 
coded temperatures relevant to the bubble phase, in units 
of e%. The model is assumed with parameters /3i — 0.07 
and Ai = 0.45 A^o, where No = p F M/ir 2 . 



is the sum of a FC part given by 

vj 

Pf(T) = — [ n (p) [1 - no(p)] dp ~ (11) 



and a noncondensate part P n (T) = P> (T) + P< (T) con- 
sisting of two terms. Defining 



P„(T;a,6) = — 

PF 



$2 



n( fl 



(12) 



Outside the FC domain, £(p) still remains flat (see 
Figs. 3, 4), and contributions to thermodynamic prop- 
erties from the regions adjacent to pi and p/ also play 
significant roles in their NFL behavior. Indeed, consider 
for example the static magnetic susceptibility 



X (T) = -Ai|lIo(w = 0)[1 - .gon (w = 0)]" 



(10) 



Here \xb denotes the Bohr magneton |2JJ and go the ze- 
roth harmonic of the spin-spin interaction function, while 
Hq(uj = 0) = J (dn(£)/d£)du is the static particle-hole 
propagator. 

The propagator n (w = 0) = -N (P f (T)+P n (T))/T, 
where Nq is the density of states of the ideal Fermi gas, 



the two terms become P,f = P n {T\ — M>£i) and Pn = 
P n (T;£f, 00). According to Eq. (J2J, at small rj one has 



^(P^P/)M> = M T )+^(p-p/) s * + 



(13) 



with U/(T) ~ T (see Eq. ©) and s = 2 and 3, respec- 
tively, for the cases pb Pf and 7J& — pf- An analogous 
formula applies for the group velocity d^/dp outside the 
FC domain close to the point p*. Based on these results, 
algebra similar to that performed in Ref. ^| leads to 

n (cj = 0) = -Not- 1 ^ + P„r 1/s ] + const , (14) 

where r = T/e F is the dimensionless temperature. 
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FIG. 2: Same as in Fig. Q but for the model @ with param- 
eters A 2 = 3 iV and f3 2 = 0.48. 



The NFL excess A%(T, p) — x(T, p) — Xfl(p) over the 
Pauli result then acquires the form 



A X (T, P ) ~p 2 B N T- 



PfV + PnT 1 '* 



l+9o(p)N T- 1 (P f r 1 +P n T 1 / s ) 



(15) 

We see that at small rj < t 1 ^ the FC plays a minor role, 
and the NFL part of x behaves as 



(16) 



By contrast, in the case rj > t 1 ^ the FC contribution to 
Eq. I|15|) is predominant, and the magnetic susceptibility 
mimics that of a gas of localized spins 



X (T) ~ (T - Gw)- 1 



(17) 



with the Weiss temperature &w — "dovip) P- 

In heavy-fermion metals, both index regimes, i.e. s = 2 
and s = 3 0, are present, while data on the magnetic 
susceptibility of 2D liquid 3 He are compatible only with 
s = 3. Fig.[5]compares results from Eq. i|15|) for 2D liquid 
3 He with experimental data of Ref.|l| at the densities p = 
0.052 A and p — 0.058 A. We have made these parameter 
choices: P n — 0.2, Pf = 1, T] = (lower curve), and 
rj = 0.04 (upper curve). The theoretical results are seen 
to be in agreement with the experimental data. 

The FC contributions to other thermodynamic prop- 
erties are found by inserting the distribution no(p) 
into the corresponding Landau formulas. In particu- 
lar, the FC entropy Sf arising at T ~ Tz is Sf = 
- J2 i n o(p) lnn (p) + (1 - n (p)) ln(l - n (p)] ~ rj. This 
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FIG. 3: Single-particle spectrum ^(p) in units of e° F at the crit- 
ical temperature Tz = 3 x 10 -3 e% (left panel), occupation 
numbers n(p) (right-top panel), and £(p)/T (right-bottom 
panel), plotted versus p/pf at four color-coded temperatures 
relevant to the phase with a fermion condensate, in units of 
The model © is assumed. 
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FIG. 4: Same as in Fig. 01 but for the model @. The single- 
particle spectrum in the left panel is shown at Tz — 3 x 
W- 4 e%. 



term does not contribute to the specific heat C(T) = 
TdS(T, p)/dT; however, it does affect the thermal expan- 
sion (3(T) ~ dS(T, p)/dp giving rise to a great en- 
hancement of the Griineisen ratio 0(T) / C(T) ^observed 
at low T in several heavy-fermion metals [3, |23J . 

Our analysis has been carried out within the Landau- 



4 



10 



- 1 — i — i i i i 1 1 1 1 — i — i i i i i | 



0,1 




T (mK) 



FIG. 5: Low-temperature magnetic susceptibility of 3 He films. 
Experimental data from Ref.[l|appear as solid and open circles 
and solid squares, while solid curves trace the predictions of 
the current theory at low T. 



Migdal quasiparticle picture where damping 7 of single- 
particle excitations is neglected. It is applicable to sys- 
tems containing a FC, provided the dimensionless damp- 
ing rate 7(T)/T remains small at T ~ Tz- To establish 
this property, we adopt the standard formula |l7| 

7 (£~T) ~£ /"|r 2 ( 9 ,^)|Imn (q,w)ImG fl ,(p-q,^-e)^ 
1 

(18) 

r being the scattering amplitude and Gr(j>,e) the re- 
tarded Green function. With a FC present, overwhelm- 
ing contributions to the collision term l|18(l are known 
to come from a region q < q c — rjpF where I mHp fq, oj) 
substantially exceeds the ideal- Fermi-gas value 20] . Fac- 
toring out |rj 2 from the integral (|18fl . a straightforward 
argument |20| yields 7(T) ~ T _1 . However, this treat- 
ment is erroneous, since it misses the suppression of the 
integral l|18[) stemming from the inequality |r(q, w)| < 
[ImIIo(<7, ui)] 1 25] . Upon inserting this inequality into 
Eq. Ijl8|l we are led to the result 



7 (T)< 



9<9c 



ImG fl (p+q,£+u;) 
|Imn (g,o;~r)| 



dqdcu 



(19) 



In the 2D electron gas this formula coincides with that 
derived in Ref. 0. Based on Eq. lj*H3|) . the putative 
behavior j(T) ~ T _1 is replaced by 



7 (T) ~ Tr, ln(l/Tj) 



(20) 



In deriving this result, we have taken account of the 
fact that in the system with a FC, the Fermi velocity 
up = pp/M* entering the corresponding expression for 
7(T) in Ref. must be replaced by the FC group ve- 
locity (d£/dp)T — T/ripp. (More details may be found 
in Ref. [24J.) Similar results are obtained in the 3D case, 
but without the logarithmic correction 111(1/77). We con- 
clude that the quasiparticle picture holds in systems with 
a FC, at least when 77 <C 1. 



Finally, we consider stability conditions for systems 
with a flat portion in the spectrum £(p). To be definite, 
we examine the spin-density wave channel and suppose 
that at the critical density rib, the associated stability 
condition, which has the form 



1 > 5o(&c)n (fc c , lj=0) = g Q (k c 



n(p)-n(p+k c ) 
£(p)-£(p+k c ) 



dv 



(21) 

is not yet violated. Beyond the density rib, the mag- 
nitude of the NFL component of LTo, proportional to 77, 
drops rapidly with increasing T (see e.g. Eq. (|T4"1) L In the 
present case this implies that violation of the inequality 
l|2ip can occur only at very low T. Furthermore, nu- 
merical calculations show that in systems with a FC, the 
NFL part of the function LTo(fc, to = 0) has a maximum at 
small nonzero k < rjpp] one expects a new ground state 
possessing a long-range magnetic superstructure, as well 
as a small ordered magn etic moment. A salient feature 
revealed in Refs. [l3j, |2ji |28j is the destruction of the flat 
portion of the spectrum £ (p) by imposition of sufficiently 
weak magnetic fields, which also kills the magnetic order- 
ing. Such a quantum phase transition has been uncovered 
in recent experiments on YbAgGe |29J . The interplay be- 
tween quantum antiferromagnetism built upon a FC, and 
the FC itself, will be the subject of a future article. 

The present investigation has revealed generic features 
of strongly interacting Fermi systems which exhibit a flat 
single-particle spectrum We have shown that the 

quasiparticle picture remains applicable in the evaluation 
of thermodynamic properties of such systems, and that 
the flattening of £(p) is a principal source of their non- 
Fermi-liquid behavior. 
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